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SUMMARY 



A closed-form thermodynamic pole approach, TPA, is developed for average description 
of the El radiative strength functions using the microcanonical ensemble for initial states. 
The method can be applied to calculate the dipole strengths in heated and cold nuclei for 
processes of the gamma-decay as well as photoabsorption in a unified way. A semiclassical 
description of the collective excitation damping in the TPA method is based on modern 
physical notion on the relaxation processes in Fermi systems. The strengths within TPA 
approach depends on the excitation energy, i.e. Brink hypothesis is violated in this method. 

The TPA calculations were compared with experimental data and calculations within the 
EGLO method which is recommended by the IAEA as the best practical model for calculation 
of the El strength. It is shown that the TPA model is able to cover a relatively wide energy 
interval, ranging from zeroth gamma-ray energy to values above GDR peak energy. It gives 
rather accurate means of simultaneous description of the 7 - decay and photoabsorption 
strength functions in the medium and heavy nuclei. The results obtained by EGLO and TPA 
approaches are almost the same at low energies e 7 ^3AfeV. In this range the EGLO and 
TPA models describe experimental data much better than the standard Lorentzian model, 
SLO, and give a non-zero temperature- dependent limit of the strength function for vanishing 
gamma-ray energy. For gamma-ray energies near neutron binding energies the calculations 
within the TPA model describe experimental data somewhat better for heavy nuclei with 
A>j 150 as compared to other closed-form approaches. The set of the best parameters for TPA 
calculations of the El strengths in heavy nuclei is determined. The further investigations 
are important to refine dependence of the collisional and fragmentation components of the 
strength function width on gamma-ray energy and mass number. 

The computer codes were created for the calculations and plotting of the radiative 
strength functions. The El strength are calculated within framework of the SLO, EGLO 
and TPA models as a function of gamma-ray energies or mass number. The codes made 
under MS-DOS and Windows 3.1X/9X operating systems are written in Fortran and Del- 
phi programming languages. An option of visual comparison between the calculations and 
experimental data is included. 
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INTRODUCTION 



Gamma-emission is one of the most universal channel of the nuclear decay, because 
it generally may attend any nuclear reaction. This process as well as absorption of the 
gamma-rays and electron - positron decay are described in the many-body systems by the 
radiative strength functions - ||). These functions are important for the study of the 
nuclear structure models, 7-decay mechanisms, deformation and fluctuation of the nuclear 
shape, energies and widths of the collective excitations ([|J - ||). Besides this fundamental 
importance from a theoretical point of view, the strength functions are necessary to generate 
the data for the energy and non- energy applications. It is critically important to have a 
simple closed-form expression for the 7 - ray strength function because in most cases this 
function is an auxiliary quantity used in calculations of different nuclear characteristics and 
processes. The theory-based approaches for 7- strength are preferred over the empirical ones 
to improve the reliability and accuracy of such calculations and to understand the physical 
sense of used parameters. 

According to Brink hypothesis flOl, |TT|1 , the Lorentzian line shape with the energy- in- 
dependent width (SLO model) is widely used for calculations of the dipole (El) radiative 
strength. This approach is most appropriate simple method for the description of the pho- 
toabsorption data on medium and heavy nuclei (|12|] - (T4J). The situation is more com- 
plicated in the case of the gamma-emission. The SLO model strongly underestimates the 



gamma-ray spectra at low energies e 7 ^, lMeV [|T|, 16]. A global description of the gamma- 
spectra by the Lorentzian can be obtained rather satisfactorily in the range 1 ^ e 7 ^ 8MeV 
but with use of the giant dipole resonance (GDR) parameters which are different from those 
based on photoabsorption data. On the whole, SLO approach overestimates the integral 
experimental data (the capture cross sections, the average radiative widths) in heavy nuclei 
( Ulif - ||2l|). The models for description of the El strengths at low energies e 7 were proposed 
in Refs. p2| 23]. An enhanced generalized Lorentzian model (EGLO) was used and analyzed 



in Refs. p4|, p5[ for a unified description of the low energy and integral data. The EGLO 
radiative strength function consists of two components (for spherical nuclei): a Lorentzian 
with the energy and temperature dependent width Tfc(e 7 , T), and finite value term from p2| 
corresponding to zero value of 7 - ray energy. An empirical expression for width Tfc(e 7 ,T) 
was used with two additional parameters. The dependence of the parameters on mass num- 
ber was obtained to fit EGLO calculations to the experimental data. Nowadays the EGLO 
method is recommended by the IAEA p5| as the best practical model for calculation of the 
dipole gamma - ray strength function when the experimental data are unavailable. 

It should be noted that the SLO and EGLO expressions for the gamma-decay strength 
function of heated nuclei are in fact the parametrizations of the experimental data, namely, 

1) these expressions are not consistent with general relation between strength function 
and the imaginary part of the response function of the heated nuclei (see [26] - j28) and 
Sect. 2); 

2) the EGLO damping width Tk has the empirical dependences on 7-energy and temper- 
ature T. It is similar to the one of the zero sound damping in the infinite fermi- liquid when 
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the collisional (two-body) dissipation is taken into account only. It is well known that an 
important contribution to the total width in heavy nuclei is given by the fragmentation (one - 
body) width which determines a redistribution of the particle- hole excitations near of the 
collective state ( f29|j - 0]). The latter component of the width is almost independent of the 



nuclear temperature. The width of the SLO model can be identified with this fragmentation 
component. 

The statistical description of the average 7 - decay strength of excited states is presented 
below using the microcanonical ensemble for initial states. The contributions to the re- 
laxation width resulting from the both interparticle collisions with retardation effects and 
fragmentation are taken into account in a semiclassical way. The dependences of the 7 - 
decay and photoabsorption strength functions on the initial excitation energy, the gamma- 
ray energy and mass number are investigated within the thermodynamic pole approximation 
(TPA method). The TPA calculations are compared with those ones within EGLO, SLO 
models and with the experimental data. It is shown that the TPA model is able to cover a 
relatively wide energy interval, ranging from zeroth gamma-ray energy to values above GDR 
peak energy. It gives rather accurate method of simultaneous description of the 7 - decay 
and photoabsorption strength functions in the medium and heavy nuclei. 

GAMMA -RAY STRENGTH FUNCTIONS IN HEATED NUCLEI 

We shall consider the radiative strength function averaged over spins of initial states 
for 7 - emission of the electric type with the energy e 7 and multipolarity A. The general 
expression for this function can be obtained from the relation for the average radiative width 



r(e 7 ) per unit of the 7 - ray energy interval | 32fl . Within statistical mechanics the width 
f (e 7 ) is defined in standard way as the quantity averaged over states with slightly different 
values of the total initial energy E and numbers of protons Z and neutrons N 

dT 

de 



rA(e 7 )= E (i) 
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where M = Y^j^iE, Z, N, J)AEAZAN is the total number of initial states; u(U, Z, N, J) 
is the density of states; U = E — E is the initial excitation energy ; Eq is the ground state 
energy. The quantities J and M are the spin of initial states and its projection on the Z 
axis, respectively; AE, AZ, AN are the small - scale intervals of the dispersion in values of 
the energy, numbers of protons and neutrons near the average values E, Z, N. 
The quantity 

dTif/dsj = dx(e 1 )B^' S(E — Ef — e 7 ) 

is the 7 - transition probability with the energy e 7 from an initial state i to the final state 
/; rf A (e 7 ) = (e 7 /ftc)( 2A+1 )87r(A + 1)/(A[(2A + l)!!] 2 ), and 

B\f = E l< J f M f E f I I JME >l 2 



M 
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is the reduced transition probability with the multipole operator Qx^ for EX radiation, 



k 

e n (A) = e(-l) x Z/A\ e p (X) = e[(A - 1) A + {-l)\Z - l)]/A\ 

The quantities e n (A) and e p (A) are the effective kinematic charge of the neutrons and protons 
in nucleus, respectively. 

Eq. can be represented in the following form 

r x (e 1 )=<D(e 1 )> /Q(U,Z,N). (2) 

Here, Vt{U,Z,N) = J^j^iU, Z, N, J) is the total density of the initial states; the symbol 
< . . . > denotes an average over the energies and numbers of the protons and neutrons with 
the unit weight functions in the intervals AE, AZ and AN, respectively; 

D(e r ) = rf A (e 7 ) E S(E - E U )5(N - N')8{Z - Z') x 

N',Z',J,M,Jf,Mf,v 

\< JfM f Ef | Q x | JME >| 2 5(E U - E f - e 7 - 7l (iV - N') - l2 {Z - Z')), (3) 

where Qx = J2^Qx^- The identical changing the arguments is made in the 5 - function 
depending on energy. The additional constants 7, defined below fix the numbers of the 
protons and neutrons. 

In the region of high excitation energies being discussed the density of states in the 
intervals AE, AN, AZ is almost constant. Therefore, we can assume that the D varies a 
little and < D >= D. In this case the quantity T\{^-y) = D/Q coincides with the width of 
the 7 - decay of states of the microcanonical ensemble with the given constants of motion 
E, Z and N. Using the integral representation of the 5 - functions 

1 r+°° 

5(x) = — — / dtexp(itx), (4) 

271 J -00 

and completeness relation for wave functions one obtains 

I r+ioo 

D{e 1 ) = / daida 2 da 3 exp(a 3 U — a\N — a 2 Z)Z({aj})T({aj}, e 7 ). (5) 

\2Tly J —ioo 

Here, Z({aj}) = Sp(exp(— j37i)) is the partition function of the grand canonical ensemble 
characterized by three constants a\, a 2 , (3 = a^, Ti. — H — jiN — j 2 Z, 7, = (Xj//3, H is the 
nuclear Hamiltonian and N, Z are the operators of the neutron and proton numbers. The 
quantity T(e 7 , {ctj}) is the mean width per unit energy of the 7 decay of the states of the 
canonical ensemble with the parameters {aij} 

r({ aj },e 7 ) = ^4^lQi,Q x ({aj},uj), u = ey/h, (6) 

7T A 
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where 

I f+oo 

%Q.(K}^) = 7T dt S P [p(n)Q x (0)Ql(t)]exp{iu;t} (7) 

/7T J — oo 

is the spectral intensity for the expectation value of the product Q\(0)Q\(t) the multipole 
operator Q\(t) in the Heisenberg representation 

Q\(t) = exv\itH/h)Q x exjp[-itn/h), Q x = Qx(0), 

i.e. the time-depended correlation function for the operator Q\(t). The canonical av- 
erage Sp[p(TC) . . .} is taken over the Gibbs ensemble with the density matrix p{7i) = 
exp[-FH\/Z({ aj }). 

Taking into account the fluctuation- dissipative relation between the spectral density of 
nonsymmetrized correlation function and response function [[3^, [34]], Eq.(§) can be written 

as 

r({« j }, e 7 ) = s A (e 7 , {aj}) exp{-phu}, (8) 

where 

s A (e 7 ,{a i }) = I Q * Qx ({aj},uj)exp{{3huj} = 

7T A 

- [rf A (e 7 )/7r](l - ex V {-f3hu})- l Imxx{u, {a,}). (9) 
The quantity x\ is the linear response function given by 

X x(u, {a,}) = Ze p (\)x?(u, {«,}) + Ne n (X) X [ n) '(w, {a,}), 

^(^{aj}) = Sp(r^Y^f)Sn k )/qUt), (10) 

where 5nk(t) is the change of the single-particle density matrix n k induced by the external 
field 

VL = ^(t)e fe (A)r A ^F v (f), q^t) = q exp[-i(a; + vq)t\, rj -> +0, go < 1 

for protons (k = p ) and neutrons (k = n). 

The integral in Eq.(^) can be evaluated by means of the saddle- point method. The 
parameters {otj} of the saddle point are found from the condition of an extremum of the 
logarithm of the integrand. They are the solutions of the equations 

dS y ({aj})/da k = 0, k = 1 -r- 3, (11) 

with 

S 7 ({aj}) = hiZ({aj}) — a>iN — a 2 Z + a 3 U + \nT({aj}, e 7 ). (12) 

Next it is assumed that the dependence of the function s\ on the saddle point parameters 
{a,} is more smooth as that of the partition function Z({atj}), namely, 

| d n lns x /d n a j |<| d n lnZ/d n a j |, n = 1,2. (13) 
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This assumption is in agreement with many investigations of the linear response function 



properties in medium and heavy heated nuclei (see, for example Refs. p7| , |28 , 35| . One 
finally obtains [J* for the average radiative width Tx(e-y)'- 

rA(e ^ - \^hc) axM) mz,N) ' (14) 

where 

a A (e 7 , Tf) = I — J s x (T f , e 7 ), s x (T f , e 7 ) = (15) 

Here, for simplicity, the designations a dependence of the functions s and \ on the parameters 
«i and «2 is not indicated in Eq. flT4]) and below. The quantity £l(Uf, Z, N) is the total density 
of the final states with the energy Uf = f/ — e 7 and the temperature Ty = /?, 

n(U f ,Z,N) =exp[5 / ({« J })]/(27r) 2 | det[[o«] | 1/2 , (16) 

where det[a^] is the determinant of the matrix with the elements 

a kl = d 2 S f /da k dai, k, I = 1 3. (17) 

The parameters a±, a<i and 0:3 = 1/(3 = l/Tf are the solutions of equations of the thermo- 
dynamical state in final nucleus, i.e. they are the solutions of the Eq. ([□]) with the entropy 
of final states 

Sf({aj}) = S 7 ({aij}) — lnr({aj}, e 7 ) = \nZ({aj}) — a±N — a^Z + a^Uf 
instead of the function S" 7 . 



The relation (|14] ) is the same one as those given by detailed balance principle j|] with 
the photoabsorption cross-section o"A(e 7 ,Ty) given by Eq.([L5|). It follows from Eq.(|14D and 
(0), @ that the rate of the 7- transitions between excited states depends mainly on thermal 
fluctuations of the multipole moments in the final states. 

General form of the a\(e^,T) coincides with that one from Refs. p6f - pSfl , |37| , p8| 



obtained by making use of the averaging over the canonical ensemble with a constant tem- 
perature T. 

The emission and absorption processes for 7-rays are generally connected with different 
radiative strengths Jl|, |14| . The gamma- decay (downward) strength function / EX determines 
the 7- emission of heated nuclei. It is associated with the average radiative widths f\(e 7 ). 
The photoexcitation (upward) strength function / EX is connected with photoabsorption 



1 ' Note that the expression for the average gamma width, with use of the microcanonical distribution 
and an assumption on a weak dependence of the width on e 7 , was first considered in |3(| . In this case the 
average radiative width is equal to r>(e 7 ) = — d\(e 7 )/ Imxxiw ,T) / -K . 
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cross-section at fixed temperature T. For the dipole transitions these functions have the 
form 

and 

7 »(^^W^ (< " n (19) 



Here, spectral function .F(e 7 , T) is given by 

N7. / 9 \ 3 

....,,,„, (20) 



£(e 7 ,r) = l/[l-exp(-e 7 /T)], (21) 

where ^H^T) = 5'p(rX) M ^u(r : )<fo^)/9a;(£) is the response function of the heated nu- 
clei to the field qu,(t)r J2u ^i^(^) an d 5n^\t) = 8n p (t) — Sn n (t) is the variation of the 
isovector single-particle density matrix. In the case of the spherical nuclei, we have 



X 



Hi 



c<j,T) = 3S'p(rYio(f)5nH)/g LJ (t) with the isovector density perturbation 5n^~\t) in- 
duced by the dipole field q w (i)rYio(r) . 

Note that the 7- decay strength function depends on temperature Tf of the final states. 
This temperature is a function of the 7- ray energy in contrast to the initial states temper- 
ature T. 

In the case of cold nuclei the radiative strength functions is also connected with the 
response function by Eqs. (|l"8l) -([2"0"1) but with factor C = 1. The scaling factor £(e 7 , T), fl2l|) , 
defines the enhancement of magnitude of the radiative strength functions in heated nuclei 
with temperature T as compared to the cold nuclei. This factor can be interpreted as average 
number of the lp-lh excited states in heated system placed in an external field with energy 

hu, 

AN lp ^ lh (u,T) = —j J de l de 2 n(e l )(l-n(e 2 ))5(e 1 -e 2 + hu)= (22) 

X r+oa 

= — / dem(ei)(l - n(ei + hw)) = 1/ [1 - exp(-hw/T)] , 
nuj Jo 

where n(e) = 1/[1 — exp((e — /x)/T)] is the Fermi distribution function for occupation of the 
single-particle states. 

In order to get a simple closed-form expression for the response function, one assums 
as usual that gamma-decay is determined by the collective motion mode which excited in 
the associated photoabsorption process. Due to this the El transitions are considered as 
corresponding to the giant dipole excitations. Next the hydrodynamic model with damp- 
ing (||39||) is applied for description of the collective motion of the neutrons against the 



protons which corresponds to the GDR in the classical picture. This approach is an exten- 
sion of the Steinwedel- Jensen (SJ) model and provides a simple description of the GDR 
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excitation simultaneously with its damping. In common with the other classical hydrody- 
namics models SJ model with damping corresponds to the semiclassical description of the 
Fermi systems by means of the Landau- Vlasov kinetic equation with truncation of the Fermi 
sphere distortion by the layers of monopole and dipole multipolarities only |4(|. Note that 
the SJ model describes volume oscillations of the transition density and these oscillations 
are almost unaffected by the dynamical distortion of the Fermi surface with multipolarities 
I > 1 f41 ]. The SJ mode plays most important role in heavy nuclei |42[ . 



We make use of the expression for the induced dipole moment within the extended SJ 



model from $14.4 of Ref.[39| and combine it with the relation for classical absorption cross- 
section and Eqs.(p~5|), (|20|). Then we get the spectral function ^ r (e 7 ,T), Eqs. (|20|) , in the 
form 

^(e 7 ,T) = 8.674- 10- 7 ^a /3 2 ^ _*_) \ /rr , , MeV~ 3 , (23) 



A 



cxp( 



where 





-e 7 /ry 

tan(z) • 



(24) 



— Tf- 



e 7 T 



c2\2 



{Ye, 



Here, e n = z n /p and f n = (NZ/A)2/(z 2 — 2) are the energy and classical oscillator strength 
of the n - resonance, respectively; z n are solutions of the equation (f(z) = (z 2 — 2) ta.nz + 2z = 
0; 



(e 7 ,T) = ±/? 



0.5e, 



1/2 



+ zr] 



e 7 + 



e; + r 2 , 



(25) 



and ao = Ane 2 Ti/{mc) = 0.305; /3o = Ro/hu; u = (4b vo i/m)(ZN / A 2 ) is the isovector sound 
velocity with the volume symmetry energy coefficient b vo i entering the semi-empirical mass 
formula; j n (z) are the spherical Bessel functions. 

The quantity V in Eqs. (24) and ([25|) is damping width of the isovector velocity v = v p — 
v n , where v p , v n are velocities of the proton and neutron fluids, respectively. It determines 
the reduced friction force vr in the equation for isovector velocity, Eq 
It can be seen that this term corresponds to the expression 



14.60) fromRef.f39 



vr = 2hJ dp(p/m)J(p,r,t)/(27fhf 



(26) 



in semiclassical picture given by Landau- Vlasov equation with a source term J(p, r, t) for 
relaxation processes [0]]. The source term is taken as the sum of two components 

P3| , [PJ, J = J c + J s : the first one, J c , is the isovector collision integral with retardation 
(memory) effects and the second one, J c , is connected with fragmantation contribution to the 
damping width (isovector one - body relaxation). The latter component is described within 
the framework of the relaxation time approximation with the relaxation time r s different 
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from infinity only in the distorted layer of the Fermi surface with dipole multipolarity. The 
expression for damping width is obtained by the use of relation ( p6|) and the expression for 



isovector collision integral from f44j . The width is the energy-dependent one and has the 
following form q 

T = r(e 7 , T) = T c (e 7 , T) + V s , r c (e 7 , T) = (m*/m)h/T c (e 7 , T), (27) 

T s = (m* / 'm)h l 't s = (m* /m)k s (e 1 )T w , T w = hv/R. (28) 

Here, r c and T s are the collisional and one-body contributions to the total width, respectively. 
The quantity m* is an effective mass of nucleon; we will use m* = m. The quantities R and 
v = 3vp/4: are the nuclear radius and average velocity of the nucleon, respectively; vp is the 
Fermi velocity. 

The quantity r c (e 7 ,T) is the collisional relaxation time for the isovector dipole distor- 
tion of the Fermi surface. It is associated with two-body collisions in the heated nucleus 
which is placed in the electric field with the frequency uj = e-y/h. For the isotropic collision 



probabilities it is given by [44] 



7~~~~jr\ = — [l + ^/2nT) 2 } . (29) 

The dependence of the relaxation time r c on the energy e 7 results from memory effects 
in the collision integral and follows Landau's prescription [(4^, [4^], [|8|, |4l| . The temperature 
dependence arises from the smeared out behavior of the equilibrium distribution function 
near the Fermi momentum in the heated nuclei. 

The quantity a in (|29|) is defined by the magnitude of the in-medium neutron-proton 
cross section a{np) near the Fermi surface 

COTl/St 

a = — — -, const = = 23.514, a in MeV, a(np) in fm 2 , (30) 

a(np) 16m* 

The magnitude of the in-media cross section a(np) is usually taken proportional to a 
value of the free space cross section aj ree {np) with a factor F, 

a(np) = F ■ a free (np), (31) 

then the parameter a can be rewritten in the form 

a = a free /F = 4.7 /F, a free = 23.514/ a free (np) = 4.7, (32) 

when the value aj ree {np) = 5 fm 2 is adopted in line with Refs. |^T| . 



2 ) Here we do not use a normalization of the damping width to the value which corresponds to the 
zero-sound magnitude for the collisional width in infinite - matter with arbitrary multipolarities of the Fermi 
sphere distortion E3, pi, M . 
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The quantities a and F determine two-body contribution B c to the damping width of 
the GDR at zero temperature and the thermal relaxation time t(T) for nuclear viscosity [^] 
in heated nuclei: 



B c = T c (e 7 = E r ,T = 0)/r r = q/a = F ■ q/a free} q = E 2 r /(T r 4ir 2 ), 

(33) 

r{T)/h = a-b-T 2 , b= {5/Q)a{pn)/a « .833a free {pn) / a free ~ 1.11, 

where the _E r and T r = T(e 7 = E r ,T = 0) are respectively the GDR energy and width 
at zero temperature; a = [cr(pp) + a{nn) + 2<r(pn)]/4 is the in-media spin-isospin averaged 
nucleon-nucleon cross section near the Fermi surface; a/ ree = 3.75/m 2 free averaged cross 
section. 

The values of the B c should be located in the interval from to 1. This condition 
determines the limiting values of a and F, 

& — C^rnin = Qi F — E max = Oij ree /q Otfree/ ^min- (3 / l) 



The isovector one-body relaxation width T s in fl28|) is taken to be similar to the 
wall formula expression T w |3(| |3l], f|2| but scaled with an energy-dependent coefficient 
fcs(e 7 )||43"[ |4], |53|, |54|. The quantal calculations within framework of a simplified RPA[^3 



show significant reduction of the one-body width in comparison with the wall value, in par- 
ticular k s w 0.1 in the range where collective phonon energy exceeds the nuclear binding 
energy and k s ~ 0.7 if collective phonon energy is negligibly small. The value k s = 0.62 was 
adopted in Ref . ||54|j . 

The spectral function T given by Eqs. (|23|) , (p4j) can be written in more convinient form 

Tie,, T) = 8.674 ■ lO" 8 - °f £ w ^f^V r\ vr ^ 

1 - exp(-e 7 /T) ^ (e 2 - e 2 ) 2 + (r(e 7 , T)e 7 ) 2 

Here, u; n = f n /f\ = (zf — 2)/(z 2 — 2); K — > oo and ay. is the peak value of the photoabsorption 
cross-section 



a r = 10.0a f n=1 /T r = 8A(NZ/A)a /T r = Q.hiz(T TRK /T r , mb, (36) 

where <Jtrk = QONZ/A is the classical sum rule in Mev ■ mb. 

The first term in expansion fl3~5|) corresponds to the excitation of the GDR with the 
energy E r = e n= i and therefore the quantities j3 , e n> i can be defined in the term of the 
GDR energy as 

Po = z n=1 /E r , e n = x n -E r , x n = z n /z 1 , z x = 2.08. (37) 
The values of the parameters q n and x n at n < 4 are given in Table 1. 
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Table 1: Values of parameters w n and x. 



n 


1 


2 


3 


4 




1 


0.070 


0.028 


0.015 


•En 


1 


2.86 


4.42 


5.97 



The imaginary part of the dipole response function xa=i associated with the Eq. (|35|) 
exhibits the resonance behaviour, in which the individual resonances have a Lorentzian 
shape with energy- dependent width. In the cold nuclei the first term of the expression (|35|) 



for the Imxx=i was obtained within the random -phase approximation ||55|| . This first term 
of (|3f|) is also in close agreement with the imaginary part of the response function of the 
heated Fermi - liquid drop on an external pressure, when approximation of the dissipative 



nuclear fluid dynamics is used for description of the system [|6|. 

The 7- decay and photoexcitation dipole strength functions, (^8|),(^9|), have the same 
temperature-dependent limiting value for vanishing gamma-ray energy and it is equal to 

7 sl (e 7 = 0, T f = T) = 7 E1 (e, = 0, T) = 

= 8.674 • l0- 8 q K ^ oo a r V r TV(e 1 = 0, T)/E$, (38) 

where q K = J2n=iQn/ x t an d Qk->oo = {^/^5)zf/(zf — 2). All value of the sum qx^oo = 
1.008656 [z\ = 2.0815) is practically contained in the first term qx=i = 1- 

TESTING OF THE CLOSED-FORM MODELS FOR CALCULATIONS OF El 

STRENGTH FUNCTIONS 

Here, the calculations of the El radiative strength functions are compared within the 
framework of the SLO, EGLO models and the approach described in Sect. II. For not very 
high excitation energies the main contribution to the spectral function T results from the 
first term of the Eq.([35|). Therefore, in what follows we use this approximation, and the 
expression ( |35"D with K — 1 is referred to as the thermodynamic pole approximation |p6 
(TPA model), T = Ttpa-, 

The SLO spectral function, T = Tslo-, h as th e Lorentzian form but with the energy inde- 
pendent width r r rather than r(e 7 ,T), 

T SL0 ^,T) = 8.674 • (40) 
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The EGLO dipole spectral function is given by |23], ^5], T = Teglo-, 



FEGLo(<h, T) = 8.674 ■ 10-Vr r 



e 7 r fc (e 7 , T) 



(^-£ r 2 ) 2 + ( e7 r fe (e 7 ,T))2 



+ 0.7 



r fc (e 7 = 0,T) 



(41) 



where the energy-dependent width Ffc(e 7 ,T) is equal to 



r fe (e 7 ,T) 



k (A) + (1 - K (A)) 



E r — e 



(27rT) ; 



(42) 



where Kq{A) is the empirical factor; eo = 4.5MeV r . For the case kq(A) = 1, the quantity 
corresponds to the expression for the collisional damping width in the infinite matter. 

The values kq(A) are mainly obtained by fitting of the low energy experimental data and 
they depend on the models used for calculations of the temperature and level density. The 
two expressions for quantity Ko(A) are recommended in 



1, A < 148, 

1 + 0.09(A - 148) 2 exp (-0.18(^4 - 148)), A > 148, 



K {A) 



when the backshifted Fermi gas model [|57| (BSFG) is taken for level densities, and 



(43) 



K {A) 



1.5, A < 145, 

1.5 + 0.131(A - 145) 2 exp (-0.154(A - 145)), A > 145, 



(44) 



for level density from [58 



Below the backshifted Fermi gas model |57j is adopted and the Eq . fl43|) for k${A) is used. 
The equations for the temperatures T, Tf have the following form 



T 



1 + 4a(U — A) 
2a ' 



(45) 



Tf 



1 + Jl + 4a{U - e 7 — A) 1 + J 1 + Aa{aT 2 - T) - Aae 



2a 2a 
where A the energy shift parameter and a the level density parameter. The values of the 
a and A are taken from the data file bei jingJbsl.dat at the rigid-body value for moment 
inertia (see, the RIPL Handbook |j25|| , Ch.5), and from the global fitting in ||59|| , namely, 



a = 0.2L4 



0.87 



MeV 



-i 



A = -6.6A~ - 32 + x ■ VIA 



-0.5 



MeV, 



(46) 



when data in Ref.[]25| are absent. Here, x = 0, 1 and 2 for odd-odd, odd-even(even-odd) and 
even-even nuclides, respectively. 
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For very small temperature T and with negative values of A the BSFG model can lead to 
negative values of the initial excitation energy U = aT 2 — T + A. In this case usual Fermi- 
gas model is used for calculation of the energy, U = aT 2 . 

The values of the GDR energy, E r , and width, r r , and the peak of the El absorption 
cross-section, a r , are considered as the temperature- independent and taken from photonu- 
clear data file beijinq-qdr.dat \ 25| (when the data exist) or from the global systematics at 



zero temperature. In the last case they equal 

E r = 31.2A~ 1/3 + 20.6yT 1/6 , T r = 0.026£ r L91 , a r = 1.2 ■ 120NZ/(AnT r ), (47) 
for spherical nuclei and 

E rA = E r /(1 + 2/3/3), £ r , 2 = E r /(1 - p/3), 

r r ,i = 0.026££i 91 , r r , 2 = 0.026££; 2 91 , (48) 

0"r,l = 0"r/3, <7 r ,2 = 2(J r /3, 

for deformed nuclei, where (3 is the quadrupole deformation parameter. All deformed nuclei 
are considered as the axially symmetric spheroids with the effective quadrupole deformation 
parameters (3. For every nucleus the quantity (3 is founded as an effective quadrupole de- 
formation parameter which gives the same value of quadrupole moment (Q) as it is in the 
case when the general expression for the Q Eq.6.53) is used allowing for the deforma- 



tions of multipolarities L = 2,4,6 with parameters /?2, I3± and (3%, respectively. The values 
of the last parameters were taken from Moller.dat file for nuclear ground-state masses and 



deformations from RIPLf^g], Ch.l, (see also |61|). The quantity of (3 is evaluated by the 
relation 

(3 = (3 2 + 0.36/3| + 0.967/3 2 /3 4 + 0.328/^ + 0.023/3 2 3 - 0.021/^ + 0.499/?|/3 4 . (49) 

For axially deformed nuclei the El strengths (|T8| ) and (|1^) are the sum of two spectral 
function T with the parameters E n i, T r< i, <x rj i and E r> 2, r rj2 , o" rj2 , respectively. 

The two approaches are used for the scaling coefficient k s {e 1 ) in Eq.(^) for one-body 
isovector width T s : 

1) the energy-independent value 

A; s (e 7 ) = A; s (e 7 = E r ) = (T r - T c (e 7 = E r , T = 0))/T w , (50) 

i.e., the quantity k s (E r ) = fc s (e 7 = E r ) is obtained from fitting of the the GDR width at zero 
temperature and it defines the one-body contribution near the GDR resonance energy; 

2) the energy- dependent value in a power approximation of the form 

, , s / KiE^ + iks^-ksiE^Ke^-E^/E^, e 7 <2£ r , 
ks ^> ~ \ k.(0), 6 7 > 2E r , 

where the quantity k s (0) = k s {e 1 = 0) defines the fragmentation width at zero energy. 
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Table 2: The optimal parameter sets at given n and the least-squares deviations per one 
degree of freedom, Xtpa> f° r TPA model. 



n 


Mo) 


B c 


Xtpa 





k B {E r ) 


0.58 


10.5 


0.5 


0.3 


0.35 


9.4 


1 


0.3 


0.46 


10.2 


3 


0.1 


0.93 


15.2 


3 


0.3 


0.46 


12.6 


3 


0.7 


0.7 


10.1 


3 


1.0 


0.93 


13.8 



The sets of optimal parameters for calculations of the TPA response function widths at 
fixed degree n of the power approximation from (|5T|) are shown in Tabl.2. The last column 
in the Table 2 contains the values of the least-squares deviations per one degree of freedom 
of the TPA calculations from experimental El strength functions data in (n/y) reaction from 
Kopecky.dat file of the RIPL-handbook( [p5| ) , 

imax. —53 2 

X 2 = E [(f Ei(Ai) - fEli^/aifEliA,)}] . 

i=i 

Here, fEl(Ai) and a[f El(Ai)] are respectively experimental average value of El strength 
and its uncertainty recommended for nucleus with mass number A{ corresponding to the i-th 
experimental value, i max = 53. The model strength functions / si (A) were calculated at 
excitation energies and gamma-ray energies equal to mean energy of El transitions, Eg El, 
from the data file, E{AA = e 7 (Aj) = EgEl(Ai). 

The values of deviations Xtpa are m m ost cases less than in the EGLO approach for 
which xeglo = 13.8. The magnitude of the least-square deviation is equal to Xslo = 50.6, 
i.e. overall fitting experimental data from Kopecky.dat file by the TPA and EGLO models 
is better than within the SLO model. 

The quantity B c , (|33"D, appears as more convenient than F in calculations of the a, 

a = q/B c , q = 0.02533 • E^/T r , (52) 

determining two-body component of the width by Eqs.(f27|), (P§D- The reason is that maximal 
value of the F is changed from nucleus to nucleus as opposite to B c , see Eqs.(|3"3"D and fl3"ip . 

It should be noted that the damping component T c given by the Eqs.(f^j), with a from 
(|32|), (|33D, is two-body relaxation width only in some specific medium where multipolarity of 
the Fermi sphere distortion does not exceed dipole multipolarity and two-body collisions are 
isotropic ones with energy- independent cross-sections a = <Jf ree - Therefore the parameters 



15 



B c , (|33"D, and F, fl3~l"|), are respectively two-body contribution to width and proportionality 
factor of the in-medium cross-section to free-space magnitude in the specific medium and 
their values account for the difference between real system and specific medium too. 

The dipole 7- decay strength functions / E1 considered as a function of mass num- 
ber are shown in Fig.l. The experimental data taken from Kopecky.dat file of the RIPL- 
handbook(f25|). Calculations were performed for nuclei from this data file (50 nuclei corre- 
sponding to (n,7) reaction) and at excitation energies and gamma-ray energies are equal to 
mean energy of El transitions. These energies are rather close to the corresponding neutron 
binding energies. The different lines connect the values calculated within framework of given 
model. Hereafter the values n = 0.5, k s (0) = 0.3, B c = 0.35 and n = 3, k s (0) = 0.7, B c = 0.7 
are used in TPA calculations. The values of index n are only indicated in the figure for 
short. 

As it can be seen from this figure, for gamma-ray energies near neutron binding energies 
the calculations within the TPA model describe experimental data in somewhat better way 
for heavy nuclei with A> 150 as compared with other approaches. The parameters n — 3, 
k s (0) = 0.7 B c = 0.7 can be recommended as more appropriate set in TPA calculations. 

In Fig. 2 the results of the calculations of the strength functions / E1 in 1M Nd with the 
initial excitations energy E which is equal to the neutron binding energy B n « 7.8 MeV are 
shown. The experimental data are taken from Ref . [|i~5[ . 

The results obtained by EGLO and TPA approaches are almost the same at low energies 
e 7 ^,3MeV '. In this range the EGLO and TPA models describe experimental data much 
better than the SLO model and give a non-zero temperature-dependent limit of the strength 
function for vanishing gamma-ray energy, see Eq.(^). The calculations by TPA and SLO 
models at the energies e^5MeV . lie more close to experimental data than within EGLO 
method. 

Figure 3 demonstrates the dependence of the 7 -decay TPA strength function on the 
initial excitation energy U. The El strength depends rather strongly on the energy U. It is 
usually named as a breakdown of Brink hypothesis . This violation of Brink hypothesis is 
growing with increasing excitation energy. The difference of the El strength function values 
calculated at different U is increased with decreasing 7 - energies and these deviations are 
more important for the 7 - transitions with energies under or of the order of the nuclear 
temperature T. 

In Figs. 4, 5 the comparison is shown between different approaches in the case the pho- 



toabsorption strength function f E1 , Eqs. (|I9D , (p9|) , (40) and (41) at different values of the 
temperature T = 0.01, 2MeV of absorbing nucleus UA Nd. The notations are the same as 
in Figs. 2. The experimental data are taken from Ref. fi~5ll . They correspond to (n, 7) reac- 
tion at e 7 = 6 — 8MeV and were obtained from photoabsorption cross-section in the range 



e 7 > 8MeV. 



The behaviour of the El strength functions calculated by the TPA method is almost in 
coincidence with SLO model in the vicinity of the GDR peak energy. It is mainly resulted 
from account of the one-body relaxation width T s , (p8[) , which is practically independent 
of the gamma-ray energy. Note that the SLO approach is probably the most appropriate 
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simple method for the estimation of the El photoabsorption strength for cold nuclei in the 
range of giant resonance peak energy. The strength function ~~f EX depends only weakly on 
the temperature if a magnitude of the T is much smaller than the gamma-ray energy. The 
form of the strength is rather sensitive to the excitation energy of absorbing nucleus at low 
energies of the 7 -rays. 

CONCLUSIONS 

A closed-form TPA approach is developed for average description of the El radiative 
strength functions. This method is not time consuming and is applicable for calculations of 
the statistical contribution to the dipole strengths for processes of the gamma-decay as well 
as photoabsorption with compound system formation. It has the following main features: 

1. The general expression between radiative strength function and imaginary part of 
the temperature response function is used. This relationship is based on microcanonical 
ensemble for initial excited states and it is in line with a detailed balance principle. 

2. The form of the temperature response function is taken within framework of the 
Steinwedel- Jensen hydrodynamic model with damping. The response function has the 
Lorentzian line shape (two for axially deformed nuclei) with width depending on 7- ray 
energy. The Landau- Vlasov kinetic approach with the monopole and dipole distortions of 
the Fermi sphere is employed to calculate the damping width which is proportional to friction 
coefficient of the isovector velocity of the relative motion of the protons over neutrons. 

3. Description of damping in the TPA method is based on modern physical understand- 
ing of the relaxation processes in Fermi systems. The contributions to the Lorentzian width 
resulting from the interparticle collisions as well as fragmentation component caused by in- 
teraction of particles with time dependent self-consistent mean field are included. A method 
of independent source of relaxation is employed to account for all contributions to width. 
The energy dependence of the collisional contribution is arisen from memory effects in the 
collision integral. 

4. The form of the El radiative strength function within framework of the TPA model 
is determined by both the Lorentzian shape of the response function with energy dependent 
width and an average number of the excited lp-lh states at given 7-ray energy. Shell 
structure and pairing correlations are included in phenomenological way by use of the level 
density parameters allowing for these effects. The TPA approach is characterized by a non- 
zero limit of the El strength for vanishing gamma-ray energy. It gives the temperature- 
dependent form of the strength, i.e. leads to a breakdown of Brink hypothesis. 

The comparison between calculations within TPA, EGLO and SLO models and exper- 
imental data showed that the TPA approach provides rather reliable method of a unified 
description of the 7 - decay and photoabsorption strength functions in a relatively wide en- 
ergy interval, ranging from zeroth gamma-ray energy to values above GDR peak energy. The 
TPA model will be useful for the prediction of the downward and upward radiative strength 
functions for cold and heated nuclei. The values n — 3, k s (0) = 0.7 and B c = 0.7 can be 
recommended as best suited set to calculations of the El strengths in medium and heavy 
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nuclei by the TPA model. It should be noted that a behaviour of the TPA strength func- 
tions is rather sensitive to the type of 7-ray energy dependence of scaling coefficient k s {e 1 ) in 
Eq . (p8|) for one-body isovector width. The phenomenological approximations (]50f) and([5l|) 
are currently used. The further investigations of the fragmentation width are necessary to 
refine the form of the scaling coefficient. 

The results obtained within this project were partly published in [|||, 5q| , namely, 

1) V.A.Plujko, Acta Phys.Pol. B30(1999)1383-1391. 
Draft version: http: //xxx.lanl.gov/abs/nucl-th/9809010; 

2) V.A.Plujko, Nucl.Phys. A649(1999)209c-213c. 
Draft version: http: //xxx.lanl.gov/abs/nucl-th/9809012. 

They also were reported on the following conferences 

1) The topical Conference on Giant Resonances, Varenna (Como Lake) Villa Monastero, 
1998, Italy. Book of abstracts. P.53, 56; 

2) Intern. Conf. Nucl. Phys. Close to the Barrier, The centennial of the discovery 
of Polonium and Radium, Warsaw University, Heavy Ion Laboratory, 30.06 - 4.07.1998, 
Warsaw, Poland. Abstracts. P.62, 63; 

3) International Conference on Nuclear Physics (49th Workshop on Nuclear Spectroscopy 
and and Nuclear Structure), 21-24 April 1999, Dubna. Russia. Book of abstracts. P.347, 348; 

4) Workshop on Collective excitations in Nuclei and other Finite Fermi systems, 21-24 
June 1999, Dubna. Russia; 

and were also presented on the RIKEN Symposium and workshop on Selected topics in 
nuclear collective excitations, March 20 - 24, 1999, Wako, Saitama, Japan. Abstracts. P.34, 35. 
As a part of the project, the modifications of the Kopecky.dat and Kopecky. readme files 



from Segment 6 of the RIPL Handbook [25] are made. The Kopecky.dat file containing 



experimental data base of El and Ml gamma-ray strength functions was converted into the 
computer readable version and passed to the IAEA. 

The computer codes were created for the calculations and plotting of the El radiative 
strength functions versus mass number and gamma-energy by the SLO, EGLO and TPA 
models under MS-DOS (see Appendix 1) and Windows 3.1X/9X (see Appendix 2 ) operating 
systems. They were passed with full description to the IAEA. The codes are written in 
Fortran and Delphi programming languages. An option of visual comparison between the 
calculations and experimental data is included. Numerical data output in the computer 
readable form is done. 
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APPENDIX 1 



GENERAL INFORMATION ON THE COMPUTER CODE fEl_vs_A.for 
AND THEIR ASSOCIATED PROGRAM FILES 



The fEl_vs_A.for is a FORTRAN code at DOS (MS Fortran 5.0 and above) 
(written by S.Ezhov and V.Plujko) for calculation of the El gamma-decay 
strength functions (f El=f _gd(Eg,U)=f (A) ) as a function of mass number A 
at fixed excitation energy U. A gamma-ray energy value is given by Eg= 
U*REDUCE. The calculations are performed for 50 nuclei f rom'Kopecky . dat ' 
file (RIPL) which contains experimental data base of radiative strength 
functions .The recommended data for (n, gamma) reaction (Key Reac.=l and 3) 
are used only. The radiative strengths are calculated within the framework 
of the TPA, EGLO and SLO models. 

HOW TO WORK WITH THE fEl_vs_A 

I. File ' fEl_vs_A.bat 'file should be run. DOS file fEl_vs_A.BAT is batch 
file for running of the DOS program files fEl_vs_A.exe, GNUPL0T.exe for 
calculation and plotting El gamma-decay strength functions. The DOS program 
file fEl_vs_A.exe is used for calculation of the El gamma- decay strength 
functions versus mass number. This program file is obtained by running at 
the DOS of the Fortran code fEl_vs_A.for (MS Fortran 5.0 and above). The 
free GNUPL0T.EXE program is used for showing results on display. It has the 
following copyright and permission (see, file GNUPLOT . gih) : 

" Copyright (C) 1986 - 1993 Thomas Williams, Colin Kelley" 

"Permission to use, copy, and distribute this software and its 
documentation for any purpose with or without fee is hereby granted, 
provided that the above copyright notice appear in all copies and 
that both that copyright notice and this permission notice appear 
in supporting documentation" . 

II. The regimes of the work should be set providing answers to the 
questions on display. There are the following settings: 
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* * 

* Specify the representation form * 

* of El strength function for output * 

* data in "RESULTS .DAT" file and * 

* the data plotting: * 

* * 

* 1 - fEl (MeV~-3) * 

* 2 - fEl/ (A~ (2/3)) (MeV"-3) * 

* 3 - fEl/(A~(8/3)Eg~2) (MeV~-5) * 

* * 



* Specify the values of the excitation * 

* energies U: * 

* * 

* 1- U = Eg from KOPECKY.DAT * 

* 2- U = neutron binding energies * 

* scaled by factor ENERCO * 

* 3- U = COMENE (identical value for all nuclei) * 

* * 



a) If set 2 



* Set scaling factor ENERCO * 



b) If set 3 

* Set a value of COMENE (in MeV) * 



20 



3) 



* * 

* Introduce value of IV = 1 - 2 determining * 

* type of the parameters used for * 

* calculations of the relaxation width. * 

* * 

* IV=1 -> Calculations at fixed two-body * 

* contribution, Be, to the GDR widths * 

* in cold nuclei. * 

* * 

* IV=2 -> Calculations at given value * 

* of the in-medium cross section, * 

* SIGMA_NP(in-medium)=SIGMA_NP(free)*F, * 

* i . e . , * 

* WIDTH_{coll,in}=WIDTH_{coll,free}*F. * 

* * 



a) If set IV=1 



* * 

* Setting of the Be value which is two-body * 

* contribution to the GDR width at zero * 

* temperature (Be is located between and 1) * 



a) If set IV=2 



* Choose value of F between * 

* and ',g8.2,'*) * 

*) It is the maximal value of F for nuclei from KOPECKY.DAT file. 
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4) 

* * 

* Setting of a gamma-ray energy value * 

* Eg in the form Eg= U*REDUCE by * 

* specifying of the value REDUCE * 

* * 

* Set value of REDUCE (from to 1) * 

* * 

5) Choice of the scaling factor form for 1-body contribution 

* * 

* Choose the form of the k_(s) which * 

* determines one-body contribution to * 

* to the damping width and set value KS: * 

* * 

* KS=1 -> k_(s) is independed of gamma-energy * 

* KS=2 -> k_(s)(Eg), power approximation * 

* with degree "n" * 

* * 

a) If set KS=1 

k_(s) (0)=k_(s) (Eg=Er) and n=0 

b) If set KS=2 

* * 

* Set of the k_(s)(0) value of the one-body * 

* contribution to width at zero energy * 

* * 

* * 

* Set of the n value of the degree for * 

* k_(s)(Eg) * 

* * 



22 



III. To see results on the screen one should press 'Enter', when message 
'Press enter' is appeared. The lines denoted as 'average.dat' correspond 
to the average values of reduced El strengths: 2.9e-9 (MeV"-3) for 
fEl/(A~(2/3)) and 4.2e-15 (MeV-5) for f El/ (A" (8/3)Eg~2) . 

To exit from the program one should press two times 'q' and 
then 'Enter' . 

The results of the calculations are located in output f ile' RESULTS. DAT' . 
The values of the least-squares deviations deviations per one freedom degree 
(chi2) of the calculations (f(A)) from experimental data also given in this 
file. 
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APPENDIX 2 



GENERAL INFORMATION ON GRSF PROGRAM 



The GRSF is the program file for Windows 3.1/9X (written by A.Mikulyak, 
scientific adviser V.Plujko) for calculation and plotting of the El 
radiative strength functions both for gamma-decay, f_gd(Eg,U), and 
for photoexcitation,f _ab(Eg,U) , processes. The radiative strengths with 
gamma-energy Eg for the compound nucleus with initial excitation energy U 
(temperature Ti) are calculated within the framework of the TPA, EGLO and 
SLO models. The radiative strength functions are calculated as functions 
either of gamma-ray energy (f(Eg)) or mass number (f(A)). In the second 
case El strengths are calculated for 50 nuclei from 'Kopecky.dat ' file 
(RIPL) which contains experimental data base of radiative strength 
functions. The recommended data for (n, gamma) reaction are used only. 



MENU COMMANDS DESCRIPTION 



FILE menu (contains common file management routines) 



SUBMENUS 

Load data file - Loading data file from the disk (Ctrl+0 shortcut) . 
Close data file - Closing allready opened data file. 
Save - Saving results on the disk (Ctrl+S shortcut) . 
Exit - Quit from GRSF. 



ACTIONS menu (contains calculation and plotting routines) 



SUBMENUS 

Set parameters - Setting calculation parameters and type of the 
calculated strength functions. 

Plot - Plotting calculations in the window (F9 shortcut) . 
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Clear - Clearing plot window (F10 shortcut) . 



Select data for plot - Selection/unselection of the models, data and 

files which can be plotted in the window. 

Calculate chi2 - Calculation of the least-squares deviations 

per one degree of freedom (chi2) of the calcu- 
lations from experimental El strength functions 
data in (n, gamma) reaction from KOPECKY.DAT file 
(for f=f(A) mode only); corresponding 
toolbar button is "chi2" (Ctrl+F9 shortcut). 



Calculate Be - Calculation of the two-body contribution to the 

GDR width at zero temperature at given F; 
corresponding toolbar button is "Be" 
(Ctrl+FlO shortcut) . 



OPTIONS menu (contains program options) 



SUBMENUS 

Set Data Libraries 

Frame Style - 
Set Plot Style - 



Setting experimental (and other) data libraries 
(all the GRSF libraries initially 
located in the LIB subdirectory ) . 

Selection style of the frame in the window. 

Setting data color, line and point style. Defaults: 
colors TPA - red, SLO - green, EGLO - navy; 
lines without points; for '*.dat ; files: 
colors - black; points without lines. 



Show Hints - 



Toggling hints on/off; default - on 



HELP menu (contains file with help) 

SUBMENUS 
Contents - 
About - 



Displayng main GRSF help file (Fl shortcut) 
Showing program information 
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TOOLBAR description 



Toolbar buttons duplicate corresponding menu items. There are 
hints for the each of the toolbar button. 



PERFORMING BASIC OPERATIONS WITH GRSF 



I. Setting parameters 

To set parameters, click on the button "Set parameters" 

( or "Actions / Set parameters" menu item. The form entitled 

"Set parameters" will be displayed. Necessary options 

on the first page should be selected. 

1) Choose type of strength functions as functions of 
gamma-ray energy f (Eg) or mass number f (A) ; 

2) Indicate type of used excitation parameters (initial 
temperature or excitation energy) ; 

3) Set type of strength (gamma-decay or photoabsorption) ; 

4) Set the values of the parameters Be or F, which define 
collisional contribution to the width; 

4) Set the values of the parameters k_(s)(0) = ksO and n = Degree, 
which define scaling factor of one-body isovector component of 
the width (the value n=0 is set if Degree < 0) . 

After selection options on the first page the "Continue" button 
should be pressed. New page will be displayed which depends on 
strength function type (f =f (A) or f=f(Eg)). In the case of f =f (Eg) 
you should enter necessary parameters in the corresponding edit 
boxes (nuclide, level density and GDR parameters) and then press 
"Ok" button. You can use parameters from data library (Ctrl+B key), 
parametrization (Ctrl+P key) or given by manually (see hints for 
shortcuts on this page) . In the case when library is used, the 
quadrupole deformation parameter 'beta' is automatically calculated 

by 

beta=3*(Er2-Erl)/(2*Erl+Er2) . 
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In the case of f (A) calculations the values of the excitation 
energies should be selected only: 

1) radiobutton U = EgEl means that excitation energies are 
equal to gamma-ray energies from KOPECKY.DAT file; 

2) radiobutton U = Bn means that excitation energies equal 
neutron binding energies . 

II. Plotting 

Plot style and frame attributes should be selected for plotting. 
For selection plot style it is necessary to click on the "Plot 
Style" button (or choose the "Options/Plot Style" menu item) . The 
color, line and point style should be selected for each data. 
Setting by default is the following: colors TPA - red, SLO - green, 
EGLO - navy; lines without points; for .dat files: colors -black 
with points without lines. To set frame attributes it is needed 
to press button "Frame Style" (or choose "Options/Frame Style" menu 
item) . Default frame attributes are the following: show legend, 
don't show X, Y grid lines; X axis range define automatically (f= 
f(Eg) -> < Eg < Egmax = U; f = f(A) -> A = 0-250). To display 
calculated data it is needed to click "TPA", "EGLO", "SLO" buttons 
(or select "Actions/Select Data for Plot" menu item) . Then you should 
click "Plot" button (or "Action/Plot" menu item or F9 key) and plot 
will be displayed. Press button "Clear plot" (or "Actions/Clear plot" 
menu item) if it's necessary to erase plot. Note, that after changing 
parameters or window resizing plot refresh itself. 

III. Working with data files 

To load data file it is necessary to click "Load data file" button (or 
set "File/load data file" menu item) and then needed data file in open 
dialogue window should be selected. To plot data files it is needed to 
click "Select data file for plot" button (or set "Actions/Select data 
for plot/Data files" menu item ) . The specified data file should be 
selected and then it will be displayed. To close data file menu item 
"File/Close data file" should be selected and the file be indicated. 

IV. Saving of the results 

To save results it is necessary to click "Save results to disk" button 
(or set "File/Save" menu item or Ctrl + S key). In the case of f=f(A) 
experimental data from KOPECKY.DAT will be saved together with the 
calculations. By default, the calculations are performed in 200 points. 
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Fig. 1. The E1 gamma-decay strength functions 
versus mass number A. 



31 




Fig. 2. The E1 gamma-decay strength functions 
of Nd 144 for U=Bn. 
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Fig. 3. The dependence of the gamma-decay strength functions 
at different excitation energies, n = 3.0. 
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Fig. 4. The photoabsorption strength functions in Nd , T = 0.01 MeV. 
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